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Abstract. This paper defines a theory of conformal parametrization 
of digital surfaces made of surfels equipped with a normal vector. The 
main idea is to locally project each surfel to the tangent plane, therefore 
deforming its aspect-ratio. It is a generalization of the theory known 
for polyhedral surfaces. The main difference is that the conformal ratios 
that appear are no longer real in general. It yields a generalization of the 
standard Laplacian on weighted graphs. 



1 Introduction 

Conformal parametrization of surfaces is a useful technique in image processing. 
The key notion is to identify the tangent plane of a surface to the field of the 
complex numbers in a consistent way. It allows to give a meaning to the notion 
of angles of two crossing paths on the surface, or equivalcntly to the notion of 
small circles around a point. A surface with such a complex structure is called 
a Riemann surface. A conformal, holomorphic or analytic function between two 
Riemann surfaces is a function that preserves angles and small circles. 

This notion has been tremendously successful in mathematics and engineer- 
ing; an aim of this paper is to define its discrete counterpart in the context of 
the surfel surfaces, defining for example discrete polynomials (see Fig. 6). It is a 
crucial notion in texture mapping for example: consider a particular animal skin 
that is known to contain small patterns of a given shape, like round disks; if this 
texture is rendered in a way that stretches these patterns into ovals, the picture 
will be wrongly interpreted by the viewer, the distortion being understood as 
conveying an information of tilt of the underlying surface. 

The technique has many other uses, like vector fields on surfaces, surface 
remeshing, surface recognition or surface interpolation. One of its main features 
is its rigidity: the Riemann mapping theorem tells you that a surface topologi- 
cally equivalent to a disc can be conformally mapped to the unit disc, uniquely 
up to the choice of three points. In this way, surfaces which are very different can 
be mapped to the same space where their features can be compared. There is 
much less freedom than in the case of harmonic mapping for example (kernel of 
the Laplacian, see (4)), which depends on many arbitrary choices, which are too 
numerous or too sensitive in many cases. This technique is surprisingly robust 
to changes in the bulk of the surface and the dependency on the boundary con- 
ditions can be relaxed as well [1], putting rigidity where the data is meaningful. 
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This technique has been widely used in the polyhedral surfaces community [1— 
7]. In this paper we will describe its adaptation to the case of digital surfaces 
made of surfels, square boundaries of voxels in Z 3 that constitute a simple com- 
binatorial surface where each edgel belongs to one or two surfels. We develop in 
this article the theory and algorithms needed for an actual implementation on 
computers of these notions to the context of surfel surfaces. 

The additional information that we use in order to give the digital surface a 
conformal structure is the data of the normal direction [8-10]. 

In the first section 3, we will present how this information encodes a non real 
discrete conformal structure. Then, in Sec. 4 we will recall elements of dc Rham 
cohomology and apply them in Sec. 5 to define a discrete Hodge star from this 
conformal structure. How this leads to a theory of discrete Riemann surfaces 
will be explained in Sec. 6-7, first recalling the real case, then generalizing to the 
complex case, which is the main technical result of this paper. 

2 Previous work 

Conformal maps are present in a lot of areas. They are called analytic, holomor- 
phic, meromorphic or monodriffic. They are related to harmonic maps because 
they are in particular complex harmonic. In the context of discrete geometry 
processing, they were used mainly by the polyhedral community, for example 
Desbrun and al. in [1, 2] or Gu and Yau [4-7]. Circle packings have been used as 
well to approximate conformal maps, see [3] and references therein. 

3 Conformal structure 

We show here how a surfel surface equipped with normals defines a discrete 
conformal structure and gives a geometric interpretation to holomorphic maps. 

A discrete object is a set of points in Z 3 , each center of its Voronoi cell called 
voxel. A voxel is a cube of unit side, its six faces are called surfels. A digital 
surface E made of surfels is a connected set of surfels. We will restrict ourselves 
to surfaces such that every edge in E belongs to at most two surfels [11]. The 
edges that belong to only one surfel are called boundary edges. Let us call (the 
indices stand for dimensions) (0o, Oi, O2) the sets of vertices, edges and surfels 
of this cellular decomposition of the surface E. 

Note that this cellular decomposition is bipartite: there exists a bicoloring 
of its vertices, that can be colored whether black or white, no adjacent vertices 
have the same color. We consider the surfels diagonals, their end points share 
the same color, forming dual black and white diagonals. 

We call r the 1-skclcton graph, whose vertices r are given by the black 
vertices and its edges I~i by the black diagonals. It can be completed into a 
cellular decomposition of the surface by adding faces J2 for each white vertex. 
Similarly we define its Poincare dual r* composed of the white vertices, white 
diagonals and faces associated with black vertices. We will refer to r as the primal 
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Fig. 1. A surfel surface 



(black) graph and to r* as its dual (white) graph. The fc-cells of r arc nothing 
else than the 2 — k cells of T*: A ~ r 2 *_ fc . 




Fig. 2. The cellular decomposition r associated with black vertices 



The data of a normal direction at each surfel is a broadly used feature of 
digital surfaces [8-10]. This normal might come from a digital scanner, or be 
computed from the digital surface itself by various means on which we won't 
elaborate. These consistent normals give an orientation to the surface. 

This normal is used to project a given surfel comprising the four vertices 
(a;, y, x' , y') to the local tangent plane. This projection deforms the square into a 
parallelogram. Its diagonals are sent to segments which are no longer orthogonal 
in general. We identify the tangent plane with the complex plane, up to the choice 
of a similitude. We call Z this local map from the cellular decomposition to the 
complex numbers. Each diagonal (x,x') and (y,y f ) is now seen as a complex 
number Z(x') — Z(x), resp. Z(y') — Z(y). 

For example we can project the standard digital plane of cubes associated 
with Pq : x+y+z = onto this (constant tangent) plane Pq and get the following 
rhombi pattern. 



4 Ch. Mercat 




Fig. 4. The digital plane x + y + z — projected. Note that F is the hexagonal lattice, 
r* its triangular dual. 



We then associate to each diagonal (x, x') G A the (possibly infinite) complex 
ratio i p of the dual diagonal by the primal diagonal, as complex numbers. 

IP ( X x') - Z{y>) Z{V) 
ip[x,x).- z{x/) _ z{x) 

This ratio clearly does not depend on the choice of identification between the 
tangent plane and the field of complex numbers. 

We will prefer the number p to the ratio i p and still call it abusively the 
ratio. This number does not depend on the orientation of the edge, for p{x, x') = 
p(x',x). The number p is real whenever the normal is orthogonal to (at least) 
one of the diagonals, that is to say, when the two diagonals are orthogonal when 
projected to the tangent plane. We will call this eventuality the real case. It is so, 
for example, in the standard plane case in Fig. 4 where its value is the constant 
p hex = tan(|) = 1/V3. The flat square grid 1? is associated with the constant 
1. Large or small values away from 1 appear whenever the surfels are flattened 
away from the square aspect ratio. The complex valucss appear when the surfel 
is slanted away from the orthogonal conformation. 
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P-\ p — 1 p — 1 p—l+i 



Fig. 5. Several surfel conformations and the associated ratio p. 

We call this data of a graph r, whose edges are equipped with a complex 
number p, a discrete conformal structure, or a discrete Riemann surface. 

We equip the dual edge (y, y') £ r* of the complex constant p(e*) — 1/ p(e). 
In the example Fig. 4, its value is the constant ptri = tan( : |) = \/3 = l//9hex- 

We define a map / : Oo ~~ * C as discrete holomorphic with respect to p if and 
only if, it respects the ratio for each surfel: 

V (*, y, x', y') G 2 , f(y') ~ f(y) = * P(x, x') (f(x') - /(*)) . 

In the continuous complex analytic theory, a holomorphic function / is a 
complex function from the complex plane to itself, which is complex differen- 
tiable; that is to say, it is recognized by the fact that its action on a neighborhood 
around a point zq £ C is locally a similitude z i— > az + b where b = /(zq) and 
a = f'(zo), the derivative of / at zq with respect to the complex variable z. It 
sends little circles to little circles. 

In the same fashion, a discrete conformal map sends each quadrilateral surfel 
to another quadrilateral whose diagonals have the same ratio p and it can be 
pictured as sending polygons that form the double A into similar polygons that 
still fit together. For example, the hexagonal and triangular lattices, hidden 
in Fig. 4, can actually be drawn on the same picture by simply joining the 
middles of the edgels together. The hexagons and triangles are simply both 
shrunk by a factor half and fit together at their vertices. A discrete conformal 
map is recognizable in the sense that it sends each of these regular hexagons 
and equilateral triangles to polygons of the same shape, touching with the same 
combinatorics. The discrete derivative is encoded in how much each polygon has 
been inflated or shrunk, and rotated. 




Fig. 6. The discrete version of the map z t—> z 3 in the hexagonal/triangular case. 
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4 de Rham Cohomology 

In the continuous theory of surfaces, the notion of complex structure relies on the 
existence of an operator on 1-forms such that * 2 = —Id, called the Hodge star. In 
orthonormal local coordinates (x, y), it is defined on 1-forms by *(/ dx + g dy) = 
—g dx + fdy. The discrete analogous of these local coordinates are given by these 
pairs of dual diagonals {x,x') and {y,y'). 

In order to follow further this analogy, we first have to define the spaces of 
discrete functions and discrete forms such as dx,dy. This is done in the theory 
of de-Rham cohomology [12]. 

In this section we recall elements of de Rham cohomology, with functions, 
boundary operator and forms, in which the notions of discrete analytic functions 
take place. The novelty is that we need to double everything to get the best out 
of the complex structure. 

We define A := J 1 U J 1 *, disjoint union of the two dual graphs, that we will 
call the double graph. Its vertices A = Oo are the same as the vertices of the 
surfel cellular decomposition 0, its edges A\ are the black and white diagonals, 
its faces A% is a set in bijection with its vertices. This doubling can look artificial 
at first sight but is in fact very useful in practice, allowing for nicer formulae 
because A is self-dual by construction. We could as well define complex functions 
as being real on the vertices To and pure imaginary on the faces Fi , but it is not 
as practical. The case is similar to the continuous where functions and 2-forms 
are essentially the same set but treated differently. 

The complex of chains C(A) = Co (A) C\{A) © C2(A) is the vector space 
(over R the field of reals) spanned by vertices, edges and faces. It is equipped 
with a boundary operator d : Ck{A) — > Ck-i(A), null on vertices and fulfilling 
d 2 = 0. The kernel ker d =: Z,(A) of the boundary operator are the closed 
chains or cycles. The image of the boundary operator are the exact chains. The 
homology of A is the space of cycles modulo exact chains. It encapsulates all 
there is to know about the topology of the surface. 

The dual spaces of forms are called cochains, denoted with upper indices 
for the dimension, are functions from chains to the field of complex numbers. 
C k {A) := Hom(Cfc(A),C). The space C°(A) of 0-forms is the linear span of 
functions of the vertices, 1-forms are functions of the oriented edges, 2-forms are 
functions of the faces. Coupling is denoted by functional and integral notation: 

the value of a 1-form a evaluated on an edge (x, x') € Ai will be denoted by J a, 
similarly, f(x) for a 0-form on a vertex, JJ F co for a 2-form on a face. The dual of 
the boundary operator is called the coboundary operator d : C k (A) — > C fe+1 (yl), 
defined by Stokes formula: 



where § denotes the circulation of a 1-form along a closed contour. A cocycle is 
a closed cochain and we note a £ Z k (A). The cohomology of A is the space of 
cocyclcs modulo the exact forms. 
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We define an exterior wedge product, denoted A, for 1-forms living either on 
edges Oi or on their diagonals Ai, as a 2-form living on surfels (>2- The formula 
for the latter is: 




o 



0- a p\ (1) 



(x,y,x',y') \( x ' x ') (VrV') (V,V') ( x > x ') 



The exterior derivative d is, as it should be, a derivation for the wedge product, 
for functions /, g and a 1-form a: 

d(fg) = fdg + gdf, d(fa) = dfAa + fda. 



5 Hodge star 

After having set the spaces where our operator is going to act, we show, in this 
section, how the discrete conformal structure allows to define a Hodge star, that 
is to say, an operator verifying * 2 = —Id on 1-forms. It breaks forms into holo- 
morphic and anti-holomorphic parts. Moreover, we show that this decomposition 
is robust to flips and that parallel planes have isomorphic decompositions. 

Definition 1. The Hodge star is defined on functions and 2-forms by 

* : C k {A) -> C 2 - k {A) 

C°(A)3f^*f: JJ *f 

C 2 (A) 3 w m +u : (*w)(x) 

and on 1-forms, on the dual edges (y,y') = (x, x')* € A\, given the (complex) 
discrete conformal structure p(x,x') = re 10 , by 

fi:'™) . = J_(~sin6 -i\ //..., <>\ 
\j»*a)- cos6\ r sin6 ) \^ yl) a ) ' 

Notice that the Hodge star is a real transformation. The fact that it is well 
defined relics on the fact that two dual diagonals are associated with inverse 
numbers: p{y,y') = l/p(x,x') = ^e~ l9 (using the former notation). It fulfills 
* 2 = Id for functions and 2-forms, and * 2 = —Id for 1-forms: 

We have, on the surfel (x, y, x' , y') € O2, 

1 /-sing ~l V _ 1 /sin 2 9 - 1 \ 
cos 2 6» V r sin0J ~ cos 2 6 \ sin 2 6 - l) ~ ~ 2 ' 

Because of this property, we can define it on the complexified forms and func- 
tions. It breaks the complex 1-forms into two eigenspaces associated to eigenval- 
ues —i and +i called respectively type (1,0) and (0, 1) forms. 
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Definition 2. A holomorphic form a £ Cq(A) is a closed type (1,0) form, that 
is to say, a is such that da = and *a = —ia. An anti-holomorphic form is a 
closed type (0, 1) form. 

A function f is holomorphic iff its exterior derivative df is as well holomor- 
phic, that is to say, on a surfel (x, y, x' , y') £ O2, fin') — fill) = i p(x, x') {f{x') — 

my 

A meromorphic form with a pole at F £ O2, * s a type (1,0) form, closed 
except on the surfel F. Its lack of closeness is called its residue at F. 

Holomorphic and meromorphic functions are tremendously important in math- 
ematics, they are behind all the keys on a calculator like polynomials, inversion, 
cosine, tangent, exponential, logarithm... This article defines the framework in 
which their discrete counterparts take place. 

An interesting feature of the theory is its robustness with respect to local 
moves. A discrete holomorphic map defined on a discrete Ricmann surface is 
mapped by a canonical isomorphism to the space of discrete holomorphic maps 
defined on another discrete Riemann surface linked to the original one by a series 
of flips. These flips arc called in the context of discrete conformal structures 
electrical moves. They come in three kinds, the third one being the flip, the 
others being irrelevant to our context. 




Fig. 7. The electrical moves. 



The third move corresponds to the flip, it is called the star-triangle transfor- 
mation. To three surfels (drawn as dotted lines in Fig. 7) arranged in a hexagon, 
whose diagonals form a triangle of conformal parameters p\, P2 and ^3, one 
associates a configuration of three other surfels whose diagonals form a three 
branched star with conformal parameters p\ (on the opposite side of pi) verify- 
ing 

PiPi = P1P2 + P2P3 + PsPi = — —. — r- ( 2 ) 

P1+P2 + P3 

The value of a holomorphic function at the center of an hexagon is overde- 
termined with respect to the six values on the hexagon. These values have to 
fulfill a compatibility condition, which are the same for both hexagons, therefore 
a holomorphic function defined on a discrete Riemann surface can be uniquely 
extended to another surface differing only by a flip [13]. 

This means in particular that a discrete holomorphic function defined on 
the standard plane in Fig. 4 can be followed through all its other parallel de- 
formations and is not sensitive to some added noise (flips deleting or inserting 
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extra voxels) provided the normal vector is unchanged with respect to the dis- 
crete plane value: the space of holomorphic functions on these parallel or noisy 
planes are in one-to-one correspondence. This theoretical robustness has yet to 
be experimentally observed in practice because the normal vectors are not inde- 
pendent and a noisy plane will have noisy normal vectors as well. This will be 
the subject of a forthcoming article. 

For real conformal structures, the formulae are simpler, we present them 
independently before generalizing them to the complex case. 

6 The Real Case 

We present in this section the formulae in the real case, for the Hodge star, the 
Laplacian, the scalar product and the different energies, Dirichlet and conformal. 

When the complex structure is defined by a real ratio p, for example in 
the hexagonal/triangular standard discrete plane case, then, for each surfel, the 
Hodge star takes the simpler form 

* : C k {A) -> C 2 - k {A) 

C°(A)3f^*f: jj *f 

C 1 (/l) 3 a *a : *a := 

J e 

C 2 {A) 3 ui i-> *w : (*w)(x) 

The endomorphism A := —d*d* — *d*dis, in this real case, the usual 
discrete Laplacian: Its formula on a function at a vertex x € Tq with neighbours 
Xi, . . . , xy £ Jo is the usual weighted averaged difference: 

V 

W)) (*) = E p& **) (/(*) - /(**)) ■ ( 4 ) 

k=l 

The space of harmonic forms is defined as its kernel. 

Together with the Hodge star, they give rise, in the compact case, to the 
usual weighted scalar product on 1-forms: 

(a,0):= ff aA*0 = (*a, */?) = (^) = | £ p(e) fa N (5) 

The £ 2 norm of the 1-form df , called the Dirichlet energy of the function /, 
is the average of the usual Dirichlet energies on each independant graph 

E D (f) := \\df\\ 2 = (df, df) = \ E p(x,x')\f(x')-f(x)\ 2 (6) 

0,x')e/ii 

_ E D (f\ r ) + E D (f\ r .) 
2 



:=f(F*), 

-- ~p(e*) f a, (3) 
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The conformal energy of a map measures its conformality defect, relating these 
two harmonic functions. A conformal map fulfills the Cauchy-Riemann equation 

*df = -idf. (7) 

Therefore a quadratic energy whose null functions are the holomorphic ones is 

Ec(f) ■= \\\df -i*d}\\ 2 . (8) 

It is related to the Dirichlct energy through the same formula as in the continuous 
case: 

E c {f) = \{df-i*df,df-i*df) 

= WfW 2 + \\\-i*df\\ 2 + Re(df, -i*df) 



= \\df\\' z +lm ff dfAdf 

= E D (f)-2A(f) (9) 

where the area of the image of the application / in the complex plane has the 
same formula 

AS) = \ JJ o dfAdf (10) 

as in the continuous case. For a face (x, y, x' , y') € O2, the algebraic area of the 
oriented quadrilateral (^f(x), f(x'), f(y),f(y'fj is given by 

df A df = 1 Im ((/(x') - f(x))(f(y')-f(y))) 
{X,V ' X ' ,y,) =-2iA(f(x),f(x'),f(y)J(y') 



When a holomorphic reference map z : Aq —> C is chosen, a holomorphic 
(resp. anti-holomorphic) 1-form df is, locally on each pair of dual diagonals, 
proportional to dz, resp. d~z, so that the decomposition of the exterior derivative 
into holomorphic and anti-holomorphic parts yields df Adf = {\df\ 2 + |<9/| 2 ) dzA 
d~z where the derivatives naturally live on faces. 

All these concepts turn into an actual machinery that can be implemented on 
computer and a full featured theory mimicking the theory of Riemann surfaces 
to a surprisingly far extent [14]. In particular in the rhombic case, the notion 
of polynomials and exponentials, differential equations, logarithm (called the 
Green function) follow through. A striking difference is the fact that constants 
are of dimension two, the constants on the graph _T being independent from the 
constants on the dual graph r* . 



7 Non real conformal structure 

In this section we give formulae for the scalar product, the Dirichlct energy and 
the Laplacian in the complex case. Nothing more is needed in order to implement 
computation of global conformal paramctrization of surfel surfaces. 



Discrete Complex 3D- Surfaces 11 



In the complex case, the dual graphs, which are independent in the real case, 
are no longer independent and are mixed together. 

The first construction that needs an adaptation is the scalar product. The 
formula (5) still defines a positive definite scalar product which is preserved by 
*, even though the last equality must be replaced by a mixed sum over the two 
dual edges: 

<«■ ffl - 1 E Rjfej, (w«)i* 1/ + ta w«» jj)- <") 

The Dirichlet energy mixes the two dual graphs as well: 



EdU) ,= m * = I y W)-fV)\ | p(e) |2 + Im { {e)) M-M \ 

(12) 

and the Laplacian no longer splits on the two independent dual graphs: 
Given x G /1 , with dual face Xq = (yi, J/2; • • • i J/v) € ^2 and neighbours 
xi, x 2 , . . . , x v G yl , with dual edges (x ,x k )* = (y k ,Vk+i) G and y y+ i = 
j/i, we have 

^(/)(*o) = E Re(p(e)) (lp(e)| 2 (/fe) " /(a)) + Im(p(e)) - f(y k ))) . 

(13) 

This Laplacian is still real and involves not only the neighbors along the 
diagonals of surfels on a weighted graph, which is the star used for the usual 
Laplacian, but all the vertices of the surfels that contains the vertex. 



8 Conclusions and acknowledgments 

This paper defines the theory of analytic functions on digital surfaces made of 
surfels. It required more than a mere adaptation of the theory known in the 
polyhedral community since the conformal parameters are in general not real. 

In forthcoming papers, this theory, already put to use in the context of poly- 
hedral surfaces, will be implemented in the context of surfel surfaces: The first 
applications that come to mind are the recognition of digital surfaces, simple 
ones and of higher topology through the computation of their discrete period 
matrices [15], the analysis of vector fields on digital surfaces allowed by the 
Hodge theorem that decomposes vector fields into rotational and divergence free 
parts, the creation of vector fields with given circulation properties, in general 
the correct discrete treatment of partial differential equations on a digital surface 
which are solved by analytic functions in the continuous case, like incompressible 
fluid dynamics for example. 

The author would like to thank the referees for useful comments, the forth- 
coming Geometric Discrete ANR project, in particular Valerie Berthc (LIRMM) 
and Remy Malgouyres (LLAIC) for discussions that brought the main idea that 
made this paper possible. 
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